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Abstract We discuss the properties of few electrons and electron-hole pairs con- 
fined in coupled semiconductor quantum dots, with emphasis on corre- 
lation effects and the role of tunneling. We shall discuss, in particular, 
exact diagonalization results for biexciton binding energy, electron-hole 
localization, magnetic-field induced Wigner molecules, and spin order- 
ing. 

1. Introduction 

The similarity between quantum dots (QDs) and natural atoms, en- 
suing from the discrete density of states, is often pointed out (Maksym 
and Chakraborty, 1990; Kastner, 1993; Ashoori, 1996; Tarucha et al., 
1996; Rontani et al., 2003). Shell structure (Ashoori, 1996; Tarucha 
et al., 1996), correlation effects (Bryant, 1987), and Kondo physics 
(Goldhaber-Gordon et al., 1998; Cronenwett et al., 1998) are among 
the most striking demonstrations. Aside from these similarities, and in 
addition to the huge technological interest in quantum systems which 
can be grown with very high control, there exist two main differences 
between natural and artificial atoms that make QDs particularly inter- 
esting also from a fundamental point of view. First, while in natural 
atoms Coulomb interactions are typically of the same order of the low- 
est energy single-particle gaps, in QDs these interactions can be made 
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larger or much larger than the latter, due to the different scaling of 
kinetic and Coulomb energy contributions with respect to confinement 
length. Secondly, also magnetic field daily attainable in laboratories are 
associated to energy scales which are comparable with single-particle 
gaps; therefore, they can be used to manipulate the quantum states and, 
in particular, the ratio between kinetic and Coulomb contributions, so 
that in QDs one can reach regimes which are unattainable in natural 
atoms. 

Coupled quantum dots, also called Artificial Molecules (AMs), extend 
to the molecular realm the similarity between quantum dots and arti- 
ficial atoms (Kouwenhoven et al., 1990; Livermore et al., 1996; Blick 
et al, 1998; Brodsky et al., 2000; Rontani et al., 2001; Pi et al., 2001). 
Inter-dot tunneling introduces a kinetic contribution which has to be 
added to the QD confinement energy, and can be tuned with respect to 
it by structure engineering and by transverse magnetic fields. In addi- 
tion, Coulomb interactions, both inter-dot and intra-dot, mix up in a 
very complicated way different electronic configurations (Slater deter- 
minants), the ratio between different contributions being a function, for 
example, of the external magnetic field strength and direction. There- 
fore, AMs constitute an interesting laboratory to study the physics of 
correlation. A specific aspect of AMs, with respect to single QDs, is that 
carrier-carrier interaction is tied with localization. Indeed, while kinetic 
energy favors derealization over e.g. two QDs forming a diatomic AM, 
Coulomb interaction favors localization in opposite dots for charges of 
the same sign; the ratio between the two contributions and, therefore, 
the localized or delocalized character of the correlated state, is controlled 
by the tunneling energy. 

From the theoretical point of view, highly correlated quantum systems 
are obviously a challenge. In many istances one is not allowed to use 
mean- field methods (Pfannkuche et al., 1993; Szafran et al., 2003). Re- 
liable results can be obtained by a Configuration Interaction approach, 
where one uses the calculated single-particle states to form a large basis 
of Slater determinants to represent the Coulomb interaction. Obviously, 
this method has limitations in the number of free-carriers which can 
be treated and of single-particle levels included in the basis set. On 
the other hand, it provides accurate results which can be quantitatively 
compared with experiments and which represent a benchmark for more 
approximate methods. 

In this work we review some theoretical results obtained for few elec- 
trons and electron- hole pairs in AMs. We specifically focus on the effects 
of carrier-carrier interactions, and the interplay of tunneling with other 
energy scales. All calculations are obtained by exact diagonalization 
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methods of the Coulomb interaction. In a typical calculations we con- 
sider N e electrons and Nh holes with the effective-mass Hamiltonian 



H = H e + H h + H eh , 



where 
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with a = e,h. Here m* and e* are the effective mass and dielectric con- 
stant, respectively, g* is the effective giromagnetic factor, [xb the Bohr 
magneton, A is the vector potential generating the magnetic field B, 
and e = |e|; all parameters are taken for the AlGaAs class of materials. 
Equations (2) neglect non-parabolicity and spin-orbit effects, but other- 
wise describe samples with realistic properties, such as layer width and 
finite band offsets, by means of the effective potentials V a (r). Typical 
samples that we shall consider consist of two identical QDs obtained 
from symmetric coupled quantum wells, grown, say, along the z direc- 
tion, of width Lw, separated by a barrier d, and with band-offset V$ 
between wells and barriers. The lateral confinement can be often taken 
as a two-dimensional (2D) parabolic potential characterized by an en- 
ergy Hujq which may assume rather different values, depending whether 
QDs are obtained by gating or by etching of a 2D electron/hole gas, by 
self-assembly, etc. The lateral confinement energy Uooq is typically of the 
order 1-10 meV and often much weaker than the confinement along the 
growth direction. 

Our numerical approach consists in mapping the single-particle terms 
H a in a real-space grid, leading to a large sparse matrix which is diago- 
nalized by Lanczos methods. Single-particle spin-orbitals thus obtained 
are then used to build a basis of Slater determinants for the iV-particle 
problem, which is then used to represent the two-body term, in the famil- 
iar Configuration Interaction approach. Coulomb matrix elements are 
calculated numerically. The ensuing matrix, which can be very large, is 
again sparse and can be diagonalized via the Lanczos method (Lehoucq 
et al., ). 
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The paper is organized as follows. In Sec. 1.2 we discuss electron- 
hole complexes, and particularly biexciton binding and localization. In 
Sec. 1.3 we discuss the phase diagram of few electrons in a magnetic 
field. In particular, in Sec. 1.3.1 we consider the effects of tunneling 
on few electrons in a large vertical field, where carriers are localized 
in the so-called Wigner molecule. Finally, in Sec. 1.3.2 we show how 
spin-ordering of two electrons can be manipulated by a magnetic field of 
arbitrary strength and direction. 

2. Electron-hole complexes in AM 

In this section we are concerned with the biexciton states in two iden- 
tical vertically-coupled QDs: the main focus is on the effects of the 
Coulomb-induced interdot correlations, which are investigated as a func- 
tion of the barrier width d (Troiani et al., 2002). The features of the 
biexcitonic states are shown to critically depend on the detailed balance 
between three phenomena: (i) the interdot tunneling, which is a single- 
particle feature, tending to delocalize the carriers and to spread them 
over the artificial molecule; (ii) the homopolar Coulomb interactions, 
which is a few-particle effect, resulting in spatial correlations between 
the two electrons (holes) that tend to minimize the repulsion energy 
by localizing the identical carriers in opposite dots; (in) the heteropolar 
Coulomb interactions, that induce spatial correlations between electrons 
and holes in order to maximize their overlap. Interestingly enough, the 
interplay between these trends can be widely tuned, either by modifying 
the structural parameters of the dots or by applying external fields, in 
order to induce non-trivial behaviours already at the few-particle level. 

In the following we fix all the physical parameters of the AM but the 
interdot distance d in order to explore the different coupling (tunneling) 
regimes. At the smallest interdot distance, the symmetric-antisymmetric 
(S-AS) splitting 2i( e,/l ) = e^'g^ — e^'^ is maximized, where es (e^s) is 
the energy of the single-particle S (AS) level. For the sample of Fig. 
1, for example, where d is varied in the range 1 3 nm, the tunnel- 
ing splitting amounts to about 12 and 6 meV for the electrons and the 
holes, respectively, at d = 1 nm. As a result, both electrons are frozen 
in the S state and the spatial distribution of each electron is uncorre- 
cted both with respect to the other electron position and to that of 
the two holes. Due to their smaller effective mass, the holes tunnel less 
efficiently than the electrons (t^ < t^): therefore the interdot spatial 
correlation, which requires the occupation of the AS states, is favored 
and acts in such a way that the carriers are always localized in different 
dots. This is shown in the left-hand insets of Fig. l(b-d), where we plot 
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Figure 1. Biexciton energies (a) and 
mean distances between carriers (b-d) 
as a function of the interdot distance 
d. The plots refer to the four states 
lowest in energy, with even (continu- 
ous lines) or odd (dotted lines) par- 
ity. The insets (b-d) show the depen- 
dence of the pair-correlation functions 
on the relative position r = ri — r-2 of 
the two carriers: in particular, we take 
y = and show the dependence on 
the z (horizontal axis) and x (vertical 
axis) coordinates. The two vertically- 
coupled GaAs/AlGaAs QDs are mod- 
eled by means of a confinement po- 
tential which is double-well like in the 
growth direction (the well width and 
band offsets are Lw = 10 nm and 



(e,h) 



v, 



= 400, 215 meV respectively) 
and parabolic in the plane (tkjg e ' = 
20,3.5 meV, m\ e:h) = 0.067,0.38 mo). 
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the dependence of the pair-correlation on the relative position between 
the carriers. By progressively increasing the barrier width, the effects 
of the homopolar Coulomb interactions between the carriers continu- 
ously increase, whereas electrons and holes are not correlated with each 
other. In this kind of configuration, resulting in a factorized wavefunc- 
tion V( r ei,r e 2,r /l i,r h2 ) ^ <^ (e) ( r ei> r e2 )<?!) ( ' l) (r w , r h2 ), the probability of 
having double occupancies (two electrons or two holes) in each dot is 
strongly suppressed, i.e., the two excitons are localized in different dots. 
In approaching the weak-coupling regime, the biexciton groundstate un- 
dergoes a rapid transition towards a maximally correlated configuration, 
where all the carriers are localized in the same dot (see the right-hand 
insets of Fig. l(b-d)). The reason why this arrangement is energeti- 
cally favored as compared to the previously discussed one is entirely 
related to the in-plane correlations between the carriers. In fact, due 
to the substantial symmetry of the electron and hole wavefunctions, the 
Coulomb energy of the two configurations is the same for both types 
of particles in the mean- field limit: the increase in the Coulomb repul- 
sion arising from the stronger degree of localization is cancelled by that 
of the Coulomb attraction. The occupation of the higher Fock-Darwin 
states, however, gives rise to additional (in-plane) spatial correlations, 
such as the ones that cause the biexciton binding energies in single QDs. 
This kind of interaction, whose closest classical analogue is an induced 
dipole-induced dipole force, is a short range one and is therefore unef- 
fective as far as the two excitons are localized in different dots. In Fig. 
1 we plot the energy levels [panel (a)] and the average distance between 
each pair of carriers [panels (b)-(d)] for the four lowest biexciton states: 
the continuous (dotted) lines correspond to the two states of even (odd) 
symmetry. Two features emerge from the plots: (i) the transition be- 
tween the weakly- and the highly-correlated configurations occurs quite 
abruptly for d ~ 2 nm; (ii) in a limited interval around this value each 
of the four states swaps its features with state of equal symmetry, and 
correspondingly an anticrossing occurs between their energy levels. 

The interdot correlation is therefore seen to play a crucial role in de- 
termining the carrier localization in artificial molecules: these exhibit a 
fully three-dimensional nature, where novel behaviors occur as compared 
to the single dots. Indeed, such features have to be taken into account 
within the design of exciton-based quantum computation schemes and 
devices (Troiani et al., 2000), where the double occupancies are known 
to spoil the required (tensorial) Hilbert-space structure. 
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Figure 2. (a) Sketch of the energy levels in a 2D parabolic potential for an AM 
at zero magnetic field. s,p,d label the value of the single-particle magnetic moment 
consistent with the corresponding atomic notation. S and AS stand for Symmetric 
and AntiSymmetric levels arising from two coupled QDs with tunneling energy 2t. 
(b) Fock-Darwin levels (Eq. 3) for an AM with L w = 10 nm, d — 3 nm, Vo = 300 
meV, and Tiuio = 10 meV. 

3. Interacting electrons in AM in a magnetic 
fields 

In this section we shall consider the effects of a magnetic field of 
arbitrary strength and direction on the few-electron states in AMs. Let 
us first summarize a few well-known effects of a field which is parallel to 
the axis of the AM. 

Figure 2(a) shows the single-particle levels at zero field, separated by 
hujQ, ensuing from the 2D parabolic potential of the in-plane confine- 
ment. The two sets of shells are separated by the tunneling energy. As 
a vertical field is switched on, the Hamiltonian can still be analytically 
solved (Jacak et al., 1998), and gives rise to the Fock-Darwin (FD) levels 

4 m = £i + hn(2n + \m\ + 1) - (huj c /2)m. (3) 

Here n is the principal quantum number describing the radial distribu- 
tion of the wavefunction, and m the azimuthal quantum number labelling 
the angular momentum which is conserved in this cylindrically symmet- 
ric configuration. The oscillator frequency £1 = (lUq + u^/4) 1 / 2 , with the 
cyclotron frequency uj c = eB/m*c, shows the competition between the 
confinement energy gaps, hujo, and the field-induced gaps, Hu c . The en- 
ergy £j is the confinement energy along the growth direction. For AMs 
with sufficiently narrow symmetric quantum well we can limit our con- 
siderations to two levels, which are the ground S and AS states (i = S, 
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Figure 3. Sketch of the ground state phase diagram for an AM with respect to field 
strength and direction, and tunneling. 



AS). The Fock-Darwin states for an AM are shown with lines in Fig. 
2(b). At zero field one recovers the s,p,... shells of Fig. 2(a). At finite 
field, the ±m orbital degeneracy splits, and at large fields the levels with 
highest m tend to form highly degenerate Landau levels. 

Figure 3 shows a qualitative sketch of the ground state phase diagram 
of few-electrons in a AM which will be discussed in the next sections. 
Let us focus, for the time being, on the vertical field effects. When a 
field is applied parallel to the growth axis, the effect is to squeeze the 
states in the QD plane. By doing so, Coulomb energy increases and, 
at sufficiently large field, the system becomes spin-polarized in order to 
gain in exchange energy. At even larger fields, kinetic energy is quenched, 
since all single-particle levels tend to become degenerate [Fig. 2(b)], and 
correlation dominates the system. Accordingly, the electron system first 
becomes a spin-polarized confined Fermi sea, laterally delocalized over 
the QDs. Then, at very high fields, the systems transforms to a confined 
Wigner crystal, in which the carriers are localized in the plane of the 
QD; the latter regime is often called Wigner molecule (Maksym et al., 
2000). 

There are two obvious directions along which the phase diagram of a 
AM can be extended: the tunneling energy and an in-plane component 
of the magnetic field, as shown in Fig. 3. These two regions will be 
discussed next. 
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Figure 4- Ground-state energy (left axis) and in-plane average radius (g) (right axis) 
vs inter-dot distance d for six electrons at B = 25 T. Sample parameters are: huJo = 
3.7 meV, Vo — 250 meV, and Lw = 12 nm. Insets show the electron arrangements in 
the different phases. 

3.1 Wigner molecules and tunneling 

As an example of the effects of tunneling in the high field regime, 
we show in Fig. 4 the calculated ground-state energy vs the inter-dot 
distance d for six electrons (Rontani et al., 2002). According to the pre- 
vious discussion, we assume that at this field carriers are spin polarized; 
furthermore, at the high field considered here, electrons in a single QD 
are expected to be localized, and quantum fluctuation to play a minor 
effect. In the inset we also show the geometrical configuration of the 
localized carriers which will be discussed in more detail in Fig. 5. At 
small distances the energy increases with d (phase I), because the ki- 
netic energy exponentially grows due to the progressive localization of 
the wavefunction into the dots: electrons occupy only S orbitals, whose 
energies increase. At these small d the AM behaves as a single QD, and 
electrons sit at the vertices of a centered pentagon, which minimizes the 
Coulomb energy, as predicted by a classical static calculation (Partoens 
et al., 1997). There is another obvious regime: when the barrier is large, 
the six electrons sit, three per QD, at the vertices of two triangles stag- 
gered by 60° (Phase III) . Close to d = 5 nm inter-dot tunneling stabilizes 
Phase II, in which electrons sit at the vertices of an hexagon. The figure 
also shows the calculated average radius (g) [g = (x,y)]; Fig. 4 shows 
that these states are incompressible, in the same sense as Laughlin's 
states of the FQHE (Laughlin, 1983). Indeed, varying d acts like an 
external pressure applied in the z direction, forcing the wavefunction to 
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change: however, due to a cusp-like structure of the energy spectrum 
(Maksym and Chakraborty, 1990), this happens only in a discontinuous 
way, except for Phase II. We show below that Phase II has very special 
properties and is stabilized by tunneling fluctuations. 




(p (rad) 

Figure 5. Angular correlation function of the three phases of Fig. 4. The three 
phases correspond to interdot distances d=2, 4.6, 8 nm, respectively. 

To analyze the ground state of the artificial molecule in the different 
regimes, we show in Fig. 5 the pair correlation function P(g, z; g , zq) = 
Hi^j (${q - Qi)K z ~ z i)$(Qo ~ Qj)$(zo ~ Zj)) /N e (N e - 1) (the average 
is on the ground state). Figure 5 shows P(g, z; g , zq) along a circle in the 
same dot (solid line) or in the opposite dot (dashed line) with respect to 
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the position of a reference electron, taken at the maximum of its charge 
density, (Qq,zq). The right column shows the electron arrangement in 
the QDs as inferred by the maxima of P(g, z; g , zq). 

At small d (Phase I) the whole system is coherent, i.e., it behaves 
as a unique QD; indeed, the correlation function peaks, forming the 
outer shell of the centered pentagon, have the same height in both QDs. 
At intermediate values of the tunneling energy (Phase II) the peaks, 
corresponding to the vertices of a regular hexagon, have different heights. 
Finally, when d is sufficiently large (Phase III) , the structure evolves into 
two isolated dots coupled only via Coulomb interaction; accordingly, the 
peaks are again of the same height, but shifted by 60° . 

It is important to note from Fig. 5 that Phase I and III are strongly 
localized phases, as demonstrated by the high peak-to-valley ratio of 
the correlation function, and quantum fluctuations play a minor role; 
therefore, electron configurations are basically determined by Coulomb 
interactions, and have completely classical counterparts (Partoens et al., 
1997). On the contrary, in Phase II tunneling fluctuations prevent elec- 
trons from localizing and the configuration has a "liquid" character. 
Such phase cannot be explained in term of Coulomb interactions solely 
and, in fact, the hexagonal arrangement shown in Fig. 5 is classically 
unstable. Therefore, tunneling fluctuations may induce melting of the 
otherwise well localized Wigner molecule in the high field regime, and 
induce reentrant liquid phases, as schematically indicated in Fig. 3. A 
discussion of the possible experimental signatures of the different phases 
in inelastic light scattering experiments can be found in (Rontani et al., 
2002). 

3.2 Effects of an in-plane magnetic field 

Next we discuss the effect of a finite in-plane component of the field, 
By. To show the effect on the single-particle states, let us suppose that 
the total field is tilted by an angle 6 with respect to the growth direction, 
so that Bh ^ 0. As shown in Fig. 6, with increasing 6 energy levels tend 
to follow the FD states backward, which are shown for comparison, since 
of course B± decreases with increasing angle; in addition, however, the 
splitting between S and AS levels decreases. This shows that an in- 
plane component of the field suppresses tunneling. Note that this effect 
is larger for higher levels. Note also that here the S/AS labelling is 
used for brevity; obviously, for a general field direction with respect to 
the tunneling direction wavefunctions do not have a well defined S/AS 
symmetry. It is important to stress that deviations from the FD states 
are expected because the energy scale associated with the in-plane field is 
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Figure 6. Evolution of single particle states in a AM when the field is rotated from 
the vertical direction by an angle 8. Sample parameters are: L w = 10 nm, d = 3 
nm, Vb = 300 meV, Tiuo = 10 meV. Solid lines represent the Fock-Darwin states £ % nm 
induced by a strictly vertical field. 




(a) (b) 

Figure 7. (a) Sketch of the energy contributions to the tunneling-induced spin-spin 
interaction, (b) Two-electron levels, with indication of the main component of the 
wavefunctions in terms of S (left boxes) and AS (right boxes) states. Parameters are: 
L w — 10 nm, d = 3 nm, Vo = 300 meV, and hu>o = 10 meV. 

comparable to the tunneling gap; for single QDs, for example, reasonable 
in-plane field will not affect the FD states, since the in-plane field energy 
scale is typically much smaller than the single-particle gaps induced by 
the quantum well confinement. 
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In AMs carriers sitting on either dot are not only electrostatically 
coupled, but also have their spin interlaced when tunneling is allowed 
(Burkard et al., 2000). This is sketched in Fig. 7(a). For two electrons 
in a singlet state it is possible to tunnel into the same dot. By doing so, 
they gain the tunneling energy t; this may compensate for the loss in the 
Coulomb energy U. This process is forbidden for two electrons in the 
triplet state by Pauli blocking. Different spin orderings, therefore, are 
associated to an exchange energy, J, which to lowest order in pertubation 
theory is J a t 2 /U. While in real molecules J is fixed by the bond 
length, in AMs it is possible to tune almost all energy scales by sample 
engineering and external fields. Since an in-plane field affects tunneling, 
we expect the exchange energy to be affected by an in-plane field as well. 

Guided by the above considerations, we next consider the two-electron 
system (?; Bellucci et al., ress). At low vertical fields the ground state 
of single and coupled QDs is known to be a singlet state (Merkt et al., 
1991; Oh et al., 1996). In the moderate field regime, therefore, the lowest 
energy levels are nearly unaffected by the rotation except for the shift 
due to the reduction of the tunneling energy, with the singlet state being 
the lowest. 

At sufficiently high vertical field a singlet-triplet transitions take place 
at a given threshold field. Since the exchange energy is proportional to 
tunneling, we expect that the threshold fields will be lowered as By 
increases (Bellucci et al., ress). This is shown in Fig. 8. The singlet 
state is stable in the low field regime. The triplet state becomes favored 
in the large field regime, whether the field is perpendicular or parallel to 
the plane of the QD; it is should be noted, however, that this happens 
by different mechanisms whether Bj_ or B\\ is large. In the former case, 
the squeezing of the wavefunction has a Coulomb energy cost which can 
only be avoided by triplet spin order. This is analogous to single QDs. 
However, while a finite Bn would not affect very much electronic states 
in single QDs, where single-particle gaps are large, in AMs the in-plane 
field affects the S/AS gap. When this vanishes, no tunneling energy is 
lost by localizing in each dot; then the spin ordering becomes irrelevant, 
and the triplet state is favored only due to Zeeman energy [Fig. 7(b)]. 

In Figs. 7(b), 8 we also show in the insets the character of the two- 
electron wavefunction of the ground state. In the low By regime, the two 
electrons occupy only the S state, either with the s symmetry with op- 
posite spin (low field) or the s and p symmetry levels with the same spin 
orientation (high field), since a large vertical field reduces the s—p gaps. 
In the large B\\ regime, on the contrary, S and AS states become degen- 
erate, and are equally occupied by the two electrons, due to Coulomb 
correlations. 
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Figure 8. Singlet-triplet phase diagram calculated for a GaAs AM with Lw = 10 
nm, d = 3 nm, Tuoq = 4 meV. The insets show the single-particle occupation in terms 
of S (left to the dashed lines) and AS orbitals (right to the dashed lines) . At low By 
only S orbitals are occupied. 
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Figure 9. Exchange energy J = Es=i — Es=o for a a GaAs AM. Same parameters 
as in Fig. 8. 

In Fig. 9 we show the exchange energy, J, defined as the difference 
between the energy of the lowest triplet and the singlet levels, as a 
function of the in-plane field at zero vertical field. This is positive (i.e., 
the singlet is the ground state) at low fields, but rapidly decreases as the 
field increases. At large fields, the exchange energy changes sign, being 
eventually dominated by Zeeman energy. 
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